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E-mail:fernando@hot. caltech.edu 

Abstract 

This  paper  shows  that  77 2  (LQG)  performance  specifications  can  be  combined  with  struc¬ 
tured  uncertainty  in  the  system,  yielding  robustness  analysis  conditions  of  the  same  nature  and 
computational  complexity  as  the  corresponding  conditions  for  7Too  performance.  These  condi¬ 
tions  are  convex  feasibility  tests  in  terms  of  Linear  Matrix  Inequalities,  and  can  be  proven  to 
be  necessary  and  sufficient  under  the  same  conditions  as  in  the  'HO0  case. 

With  these  results,  the  tools  of  robust  control  can  be  viewed  as  coming  full  circle  to  treat 
the  problem  where  it  all  began:  guaranteeing  margins  for  LQG  regulators. 


1  Introduction 

The  advent  of  modern  control  in  the  60s  brought  a  substantial  transformation  in  control  theory, 
with  state-space  tools  and  optimal  control  offering  the  promise  of  tractable,  systematic  methods  for 
multivariable  control  design.  This  era  was  epitomized  by  the  solution  of  the  LQG  control  problem 
(see,  for  example,  [1])  ,  which  provides  an  elegant,  easily  computable  method  for  a  well-motivated 
multivariable  control  design  problem:  optimizing  the  rejection  of  white  noise  disturbances  for  a 
closed  loop  system.  It  became  increasingly  clear  in  the  late  70s  that  modern  control  unfortunately 
provided  limited  tools  to  further  treat  model  uncertainty,  a  fundamental  requirement  for  a  practical 
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feedback  theory  and  an  issue  which  was  often  better  addressed  by  the  otherwise  more  primitive 
frequency  domain  techniques  of  classical  control  [20]. 

While  LQ  state  feedback  was  shown  to  provide  stability  margin  guarantees  [37],  further  research 
led  to  a  counterexample  showing  that  full  LQG  controllers  had  none  [10].  This  motivated  efforts 
to  reconcile  LQG  with  classical  methods  [11],  with  some  initial  success  in  providing  a  robust  LQG- 
based  methodology  [9].  The  most  popular  development  was  LQG/LTR  [12,  2,  40],  a  multivariable 
version  of  classical  loopshaping  using  LQG  machinery.  The  problem  of  adding  plant  uncertainty 
directly  to  LQG  remained  unsolved,  however,  and  ultimately  these  efforts  pointed  in  other  directions 
[12],  particularly  toward  (structured)  singular  values  and  related  methods  [14,  41]. 

At  about  the  same  time  as  the  critique  of  LQG  robustness  was  becoming  widely  accepted,  the 
new  performance  paradigm  of  Hoo  was  being  put  forth  [44].  It  had  close  ties  to  the  frequency 
domain  and  allowed  singular  value  robustness  conditions  to  be  treated  directly.  More  importantly, 
it  allowed  for  the  first  time  a  very  natural  and  relatively  transparent  analysis  of  robust  performance 
[13,  26].  While  Hoo  soon  replaced  LQG  (now  referred  to  as  Hi)  as  the  centerpiece  of  modern  control, 
and  research  on  Hoo  flourished  in  the  1980s,  several  developments  helped  bring  Hi  back  into  the 
picture. 

The  main  weakness  of  Hoo  is  that  modeling  signals  as  weighted  norm  balls  ignores  important 
structure,  typically  expressed  in  terms  of  spectral  or  correlation  properties,  which  are  often  features 
of  more  realistic  models  of  physical  disturbances.  Ignoring  this  structure  makes  a  worst-case  mea¬ 
sure  like  Hoo  substantially  conservative,  in  much  the  same  way  as  what  happens  when  uncertainty 
structure  is  ignored  in  singular  value  robustness  conditions;  the  recognition  of  these  limitations  led 
to  a  resurgence  of  interest  in  772  as  a  performance  measure.  The  desirable  design  specification, 
from  both  the  performance  and  uncertainty  points  of  view,  appears  to  be  in  most  cases  Robust  Hi 
performance:  rejection  of  white  signals  in  the  worst-case  over  a  set  of  plants.  (The  C\  theory  [7]  is 
another  attractive  alternative  to  Hoo  but  still  suffers  from  a  pessimistic  signal  description). 

Renewed  interest  in  Hi  performance  was  also  stimulated  by  the  striking  fact  that  the  most 
powerful  computational  solutions  for  the  Hoo  control  problem  [18,  15]  relied  on  the  same  state-space 
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tools  as  LQG.  This  led  to  a  new  research  direction  in  mixed  772/77oo  control  (see  e.g.  [5,  22,  45, 
16,  34]),  and  to  various  upper  bounds  for  the  7f2  cost  over  a  set  of  plants  (e.g.,  [42,  33,  45,  19,  17]). 

In  spite  of  these  developments,  the  robust  TL2  problem  lagged  substantially  behind  Hoo  (or 
£1),  where  a  sophisticated  set  of  tools  is  available  for  the  analysis  of  robust  performance  under 
structured  uncertainty  (see  e.g.  [26,  43,  4,  7]),  including  several  results  that  exactly  analyze  robust 
performance  with  structured  uncertainty  in  terms  of  computationally  attractive  convex  conditions 
([21,  38,  25,  35]).  No  such  results  have  previously  been  available  for  robust  H2  performance. 

This  paper  provides  the  final  step  in  the  return  of  the  H2  performance  paradigm,  casting  it 
on  an  equal  footing  with  Hoo  •  We  present  a  convex  condition  for  robust  H 2  performance  analysis 
under  structured  uncertainty,  of  a  very  similar  nature  to  the  the  corresponding  condition  for  robust 
Hoo  performance,  and  with  analogous  properties.  Computationally,  it  reduces  to  a  Linear  Matrix 
Inequality  (LMI,  see  [6])  over  frequency  which  can  be  handled  with  analogous  tools  as  in  the 
Hoo  case.  From  a  theoretical  point  of  view,  the  condition  is  shown  to  be  necessary  and  sufficient 
under  the  same  assumptions  for  the  uncertainty  as  in  the  corresponding  Hoo  conditions.  The  tools 
involved  in  proving  the  necessity  results  build  on  recent  work  in  the  Integral  Quadratic  Constraint 
(IQC)  formulation  which  has  been  mainly  applied  to  describe  uncertainty  [25],  but  can  also  be  used 
[24,  27]  for  signal  characterization.  In  this  paper  we  extend  these  methods  to  set  characterizations 
of  white  signals  based  on  statistical  tests  on  the  cumulative  spectrum  [28],  and  rely  on  infinite 
dimensional  convex  analysis  methods  to  derive  the  necessary  conditions. 

The  paper  is  organized  as  follows.  In  Section  2  some  material  on  standard  robust  control  with 
Hoo  performance  measure  is  reviewed.  Section  3  discusses  H2  norms  and  set  characterizations  of 
white  noise  signals.  In  Section  4,  the  condition  for  robust  H 2  performance  is  presented  and  proven 
to  be  necessary  and  sufficient  under  various  uncertainty  assumptions.  In  Section  5  we  remark  on 
the  computational  properties  of  this  test.  Section  6  compares  these  results  to  the  previous  work  in 
the  mixed  H2/H00  literature.  Some  remarks  on  robust  H2  synthesis  are  given  in  Section  7,  and  the 
conclusions  are  presented  in  Section  8.  Some  proofs  are  covered  in  the  Appendix.  A  partial  version 
of  these  results  was  presented  in  [31]. 
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2  Background  and  Notation 

The  results  in  this  paper  will  be  presented  for  discrete  time  systems.  Analogous  conditions  hold 
for  the  continuous  time  case,  the  details  of  which  will  be  reported  elsewhere. 


2.1  Uncertain  Systems  in  LFT  Form 


A  standard  setup  for  robustness  analysis  is  depicted  in  Figure  1,  consisting  of  a  nominal  map  M 
and  a  perturbation  A  which  enters  the  system  in  feedback  fashion;  the  overall  uncertain  system 
will  be  denoted  by  (M,  A). 


q 
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Mu  M12 

M — 

— 

M21  M22 
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Figure  1:  Uncertain  system 


M  will  be  assumed  to  be  a  finite  dimensional,  linear  time  invariant  (LTI),  stable  system.  Its 
transfer  function  is  denoted  by  M(e The  uncertainty  A  is  assumed  to  have  spatial  structure  of 
block  diagonal  form 

A  diag  [8\Iri ,  ■  - . ,  8 ],Ir ^ ,  A/,-|_i ,  • . . ,  A l+f]  (1) 

The  blocks  in  A  can  in  general  represent  real  parameters  or  dynamic  perturbations;  in  each  case, 
there  is  a  restricted  class  A  of  allowed  perturbations,  which  are  usually  assumed  normalized  to  the 
ball  of  uncertainty  Ba  =  {A  £  A  :  ||A||  <  1}  in  some  operator  norm. 

We  will  consider  I2  signal  norms:  IV;  denotes  the  space  of  squ  are- summ  able,  Cn-valued  sequences 
over  the  integers  (or  positive  integers).  These  are  identified  via  the  Fourier  transform  with  square 
integrable  functions  on  the  unit  circle  T,  with  respect  to  the  normalized  Lebesgue  measure  The 
vector  dimension  will  be  omitted  when  clear  from  context. 
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Cc(h)  denotes  the  set  of  causal,  linear,  bounded  operators  in  / 2 •  The  largest  class  of  uncertainty 
considered  here  is  the  set  of  structured  linear  time- varying  (LTV)  perturbations 

Altv  =  {A  G  £c(h)  :  A  =  dia,g[8iIri,...,8LIrL,AL+i,...,AL+F]}  (2) 

The  uncertainty  can  also  be  restricted  to  be  linear  time-invariant,  which  means  it  commutes  with 
the  unit  delay  operator  A.  This  gives  the  structured  set  ALTI  =  {A  E  ALTV  :  A  A  =  AA}. 
Some  recent  work  [35]  has  shown  it  is  useful  to  introduce  the  mildly  larger  class  of  slowly  varying 
operators,  by  defining  for  v  >  0  the  class 

A"  =  {A  £  ALTV  :  ||  A  A  -  AA||  <  v}  (3) 

of  operators  with  “variation  slower  than  vn .  For  v  =  0  we  recover  ALTI,  but  some  of  the  necessary 
conditions  will  be  proven  for  an  arbitrarily  small  v  >  0.  The  unit  balls  of  uncertainty  for  each  class 
are  denoted,  correspondingly,  Baltv,  Balti,  and  Ba^ 

The  system  of  Figure  1  is  said  to  be  robustly  stable  if  M  is  stable,  and  if  I  -  AMu  has  an 
inverse  in  Cc(h)  for  every  A  G  BA.  When  this  holds,  the  closed  loop  map  from  u  to  y  is  well 
defined  for  all  A  E  BA  and  given  by  the  Linear  Fractional  Transformation  (LFT) 

A-k  M  :=  M22  +  M2\A(I  -  M11A)~1Mi2  (4) 

A  performance  specification  can  then  be  imposed  on  the  map  A -k  M .  In  our  case  of  I2  signal 
norms,  the  standard  choice  is  the  /2-induced  norm  (which  we  call  the  Hoo  norm,  although  this  is 
an  abuse  of  notation  for  non-LTI  systems).  The  system  is  said  to  have  robust  'Hoo  performance  if 
it  is  robustly  stable,  and 

sup  ||  A  *  M||;2^/2  <  1  (5) 

A£Ba 
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2.2  Robust  Hoo  Performance  Tests 

The  main  method  for  obtaining  tractable  robust  Hoo  performance  tests  is  to  add  scalings  to  a  small 
gain  condition.  For  this  purpose  we  introduce  scaling  matrices  of  the  form 

X  —  diag  [-ST i ,  .  •  •  ,  Xjj,  X L-\-\Im\  •>••••>  XL-\-Flmjp\  (b) 


which  commute  with  the  elements  in  A.  We  will  denote  by  X  the  set  of  positive  definite,  contin¬ 
uous  scaling  functions  X(u;)  with  the  structure  (6).  The  tests  for  robust  H^  performance  can  be 
summarized  as  follows: 


Condition  1  There  exists  a  function  X{ui )  £  X  such  that  for  all  lo  £  [0,27t], 


X{u)  0 
0  I 


X(lo)  0 
0  I 


<  0 


(7) 


Since  M  is  finite  dimensional,  M(eJLU)  is  continuous  and  the  continuity  assumption  of  X(u) 
entails  no  loss  of  generality.  For  this  case  of  frequency  dependent  scales,  we  can  state  the  following: 

Proposition  1  If  Condition  1  holds  for  a  function  X{u)  E  X ,  then  the  system  [M,  A)  has  robust 
Tioo  performance  for  A  6  Balti  . 

The  previous  result  follows  by  showing  that  this  condition  provides  a  bound  for  the  structured 
singular  value  p  [14,  26],  which  is  the  exact  robustness  test  for  LTI  uncertainty.  Although  Con¬ 
dition  1  is  in  general  conservative  for  this  case,  it  remains  as  an  attractive  condition  since  exact 
computation  of  p  is  not  tractable.  Also,  computational  experience  with  bounds  such  as  those  used 
in  [4]  shows  evidence  that  the  two  tests  are  not  far  apart,  at  least  for  full  block  structures. 

Another  argument  to  support  the  claim  that  Condition  1  has  small  conservatism  is  given  by 
the  following  result  from  Poolla  and  Tikku  [35] : 
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Proposition  2  There  exists  u  >  0  such  that  the  system  (M,  A)  has  robust  Tioo  performance  for 
A  6  Ba"  if  and  only  if  there  exists  a  function  X (u>)  E  X  satisfying  Condition  1. 

Consequently,  if  one  is  willing  to  include  an  arbitrarily  small  amount  of  time  variation,  in  the 
sense  of  (3),  in  the  uncertainty,  Condition  1  characterizes  exactly  the  robust  performance  problem. 

Finally,  if  we  allow  an  unrestricted  time  variation,  robustness  analysis  is  obtained  from  Condi¬ 
tion  1  by  imposing  X(u)  to  be  a  constant  function.  The  following  result  was  shown  independently 
by  Shamma  [38]  and  Megretski  [25]: 

Proposition  3  The  system  ( M ,  A)  has  robust  TCoo  performance  with  A  E  Baltv  if  and  only  if 
there  exists  a  constant  matrix  X(w)  =  X  E  X  satisfying  Condition  1. 

The  tests  provided  by  Condition  1  amount  to  an  infinite  dimensional  convex  feasibility  problem, 
in  terms  of  a  parametrized  (by  frequency)  family  of  Linear  Matrix  Inequalities  (LMIs).  Conditions 
of  this  type  allow  for  tractable  computation,  as  will  be  discussed  in  Section  5. 


2.3  Mathematical  Preliminaries 


The  following  mathematical  facts  are  collected  here  for  ease  of  reference  (see,  e.g.,  [23]  and  [39]). 

First,  we  introduce  the  space  BV[a,  b]  of  real-valued  functions  of  bounded  variation  in  the 
interval  [a,b]  C  M.  A  function  \P(f)  is  of  bounded  variation  if 


N 

TV(fB)  :=  sup  £  | *&)  -  ¥(*,•_!)!  <  oc 

i=i 


(8) 


where  the  supremum  is  taken  over  partitions  of  [a,  b].  TF('F)  is  called  the  total  variation  of  IP. 

We  will  also  use  the  space  Cr[o,  b]  of  continuous,  real- valued  functions  on  [a,  b\,  with  the  norm 
IblU  ■=  snpte[aM\g(t)\. 

Given  *P  E  BV[a,b }  and  g  E  Cr[o,  6],  we  introduce  the  Stieltjes  integral  (see  [39]) 


/ 


g(t)d9(t) 


(9) 
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Since 


f  g(t)dV{t)  <|MU TV{V)  (10) 

J  a 

the  map  F<j  :  Cr[(i,  given  by  F,y(<?)  =  jj*  g(t)d'S(t)  defines  a  bounded  linear  functional  on 

€/;,[«,  b].  In  fact,  the  Riesz  representation  theorem  states  that  every  functional  in  the  dual  space 
Cr[o,  6]*  is  of  this  form. 

We  will  also  use  the  formula  of  integration  by  parts  for  the  Stieltjes  integral, 

[  g(t)dV(t)  =  g(b)V(b)-g(a)9(a)-  I  9(t)dg(t)  (11) 

J  a  J  a 

which  holds,  for  example,  for  $  £  BV[a,b],  g  £  Cr[o,  b}.  Furthermore,  if  g  has  an  integrable 
derivative  the  integral  on  the  right  can  be  written  as  J ^  '&(t)g'(t)dt. 

Finally,  a  key  element  in  the  proofs  of  this  paper  is  the  following  geometric  version  of  the 
Hahn-Banach  theorem,  taken  from  [23]: 

Theorem  4  Let  K, K-2  be  convex  sets  in  a  real  normed  space  V,  such  that  1C 2  has  non-empty 
interior,  and  K\  contains  no  interior  points  of  K-2-  Then  there  exists  a  bounded  functional  F  £  V* , 
r  ^  0,  and  a  real  number  a  such  that 

r(*i)  <  a  <  T{k2),  for  all  k\  £  IC\,k2  £  K, 2  (12) 


3  White  Signals  and  H2  Norms 

As  argued  in  the  introduction,  Ttoo  performance  takes  a  conservative  view  of  disturbances;  in 
many  situations  a  more  useful  performance  measure  is  given  by  the  TL2  norm  of  a  system,  which 
characterizes  the  response  to  white  signals.  For  an  LTI  system  H(e JW),  this  norm  is  defined  by 

(r  27r  duj  \  2 

Jq  trac  (13) 
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White  signals  typically  arise  in  two  situations.  One  is  as  chaotic,  high  dimensional  fluctuations 
known  as  white  noise,  which  are  usually  modeled  as  a  stationary,  uncorrelated  random  process  of 
unit  covariance  matrix;  if  such  a  signal  is  input  to  an  LTI  system,  the  variance  of  the  output  is 
given  by  \\H |||.  Another  source  of  white  signals  are  impulsive  disturbances,  or  impulsive  signals 
used  to  test  the  response  of  a  system  to  fixed  reference  signals;  the  output  energy  for  an  impulsive 
(scalar,  or  vector-valued  with  random  direction)  input  is  ||if||2-  For  more  motivation  see  [46]. 

The  objective  of  this  paper  is  to  analyze  the  effect  of  white  signals,  in  the  worst-case  for  the 
uncertain  system  given  in  Figure  1.  For  the  case  of  LTI  uncertainty,  the  closed  loop  A  *  M  is  LTI 
and  we  will  simply  analyze  for  the  worst-case  Ti-2  norm  as  defined  in  (13). 

We  will  also  want  to  consider,  however,  the  classes  ALTV  and  A "  involving  time- varying  uncer¬ 
tainty,  which  come  in  naturally  to  characterize  the  necessity  of  the  robust  performance  conditions, 
as  was  shown  in  Section  2  for  the  Hoo  case.  For  this  purpose,  it  will  be  convenient  to  describe  white 
signals  in  terms  of  a  set,  rather  than  a  random  process,  which  will  allow  the  natural  formulation  of 
worst-case  analysis  problems  over  the  uncertainty  and  over  this  set. 

A  non-stochastic  treatment  of  H 2-performance  was  in  fact  given  in  Zhou  et  al.  [45],  where  the 
classes  of  bounded  power  and  bounded  spectrum  signals  are  employed  to  motivate  both  the  7f2  and 
the  Hoo  norms.  This  formulation  is  conceptually  appealing,  but  poses  a  number  of  mathematical 
difficulties.  First  of  all,  as  noted  in  [45],  the  formalization  of  such  classes  woidd  require  limiting 
arguments  which  raise  a  number  of  technical  issues.  More  importantly,  these  classes  do  not  have  a 
rich  mathematical  structure,  which  greatly  restricts  the  applicability  of  functional  analytic  tools. 

As  a  counterpart,  the  class  of  bounded  energy  ( I2 )  signals  offers  the  rich  mathematical  structure 
of  a  Hilbert  space.  This  structure  plays  a  key  role  in  the  most  powerful  results  on  control  with 
an  TLoq  performance  measure,  and  is  equally  satisfactory  from  a  conceptual  point  of  view,  since 
bounded  power  and  bounded  energy  signals  differ  only  in  their  asymptotic  behavior.  For  this 
reason  our  treatment  of  white  signals  sets  will  be  based  on  l2-space;  this  paper  will  consider  the 
discrete  time  version,  which  is  more  straightforward.  The  same  methodology  can  be  applied  to 
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the  continuous  time  case,  as  will  be  reported  elsewhere,  and  could  also  be  used  to  provide  a  more 
complete  foundation  to  the  material  in  [45]. 

We  begin  with  the  case  of  scalar  signals  in  h-  Ideally,  a  white  signal  has  a  flat  spectrum,  i.e. 
lies  in  the  set 

W0  =  {/  G  h  ■  l/M!2  =  \\f\\l,  w  a.e.  in  [0,  2tt]}  (14) 

A  key  technical  requirement  for  the  results  of  this  paper  is  to  introduce  a  set  of  signals  which  are 
approximately  (up  to  a  small  accuracy)  white.  For  this  purpose,  we  take  the  standpoint  (developed 
extensively  in  [28])  that  such  a  notion  should  be  based  on  standard  statistical  tests  for  stochastic 
white  noise. 

More  specifically,  if  one  is  given  a  time  series  /o, . . .  ,/jv- i,  deciding  whether  it  is  a  sample  of 
white  noise  is  usually  done  based  on  the  values  of  a  chosen  statistic;  one  such  choice  is  the  sample 
autocorrelation,  and  leads  to  a  definition  for  white  noise  sets  which  was  exploited  in  [28,  29]  for 
robust  H2  analysis.  This  paper  is  based  on  a  frequency  domain  definition  for  white  noise,  which 
corresponds  to  the  so-called  Bartlett  cumulative  periodogram  test  for  time  series.  This  test  consists 
of  accumulating  the  periodogram  (squared  magnitude  of  the  Discrete  Fourier  Transform  of  the  series 
fo, . . . ,  //v— 1),  and  comparing  the  result  uniformly  with  a  linear  function.  For  more  details  see  [28]. 
Inspired  by  this,  we  consider  here  the  difference 

C/W  <15) 

between  the  cumulative  spectrum  and  a  linear  function,  and  bound  it  in  a  uniform  sense.  Define 
the  set  of  “white  up  to  accuracy  rf  signals 

wv  :=  {f  eh-  sup  \Ff(s)\  <  7]}  (16) 

•s€[0,27t] 

The  gain  of  an  LTI  single  input  system  H(e Ja;)  under  signals  in  Wo  is  easily  seen  to  be  ||ff||2. 
We  now  consider  the  worst-case  gain  under  signals  in  Wv, 

\\H\\Wv  :=  sup{||tf /||2  :feWv,  ||/||2  <  1}  (17) 
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Lemma  5  Let  Y(u)  E  BV[ 0,  2tt]  .  If  f  E  Wv,  then 


f 


2' Y(u)\f(u)\^ 


<  T]  TV(Y) 


(18) 


Proof:  Defining  Fj(s )  as  in  (15)  (Ff( 0)  =  Fj(2n)  —  0),  an  integration  by  parts  yields 

r  j.  .  r  27T 

Jo  y(w)(|/(a;)|2-||/||2)-  =  -yo  Ff(oj)dY(u)  (19) 

Since  /  E  14^,  then  ||i'/(a;)||00  =  supw  |Fy(w)|  <  rj,  so  (10)  implies  that  the  right  hand  side  of  (19) 
can  be  bounded  by  r]  TV(Y). 

a 


A  consequence  of  this  Lemma  (picking  Y(u)  =  \H(e^)\2)  is  that  for  an  LTI  system  H, 

\\H\\l  <  \\H\\2Wri  <  \\H\\l  +  TV(\H\2)  (20) 

which  implies  that  under  the  mild  assumption  that  \H(e^u)\2  E  BV[0,2ir], 

\\H\\wv  ^  \\H\\2  (21) 


We  now  analyze  the  multivariable  case;  vector- valued  white  noise  is  characterized  by  having  a 
matrix  spectrum  equal  to  a  constant  times  the  identity  matrix  across  frequency.  In  the  I2  setting, 
ideally  white  signals  do  not  appear  since  the  spectrum  is  always  a  rank  one  matrix,  but 

for  r]  >  0  we  can  define  the  set  of  approximately  white  signals  in  Iff 


W? 


f: 


/(^)/(w) 


dw 


<  V 


(22) 


2ir  2tt  m 

where  the  infinity  norm  of  a  continuous  matrix  function  is  taken  to  be  the  maximum  across  the 
coordinates  of  the  supremum  norm.  It  is  easy  to  show  that  for  every  77  >  0,  W™  is  non-trivial. 
With  this  definition,  an  extension  of  Lemma  5  can  be  written,  which  leads  for  LTI  systems  to 


ll-S'llw”*  sup{||if/||  :  /  E  W™ ,  ~\\j  M2 


<  1}  \\H ||2 


(23) 
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The  normalization  by  a  factor  in  the  input  norm  is  done  for  convenience,  since  if  )* 

approximates  the  ideal  unit  white  spectrum  Im ,  ||/||2  is  approximately  yjm. 

If  one  is  interested  in  analyzing  white  noise  rejection  for  systems  which  are  not  LTI,  the  H2 
norm  (13)  is  no  longer  meaningful,  but  it  is  natural  give  a  definition  based  on  ||Ifj|wm  as  in  (23): 

prill  -  lim  l|tf|k„»  (24) 

v— '•0+  ’ 

This  system  measure  (a  seminorm)  captures  the  response  to  signals  of  flat  spectrum,  the  interesting 
object  from  the  point  of  view  of  applications,  and  extends  the  LTI  definition. 


4  Main  Results 


In  this  section  we  provide  conditions  for  Robust  H2  performance  for  the  uncertain  system  of  Figure 
1,  which  are  analogous  to  those  presented  in  Section  2.2  for  Robust  Hoo  performance.  We  now 
state  an  analysis  test,  which  is  a  convex  feasibility  condition  on  the  unknowns  X ,  Y 


Condition  2  There  exists  X(u)  £  X,  and  a  matrix  function  Y(u)  =  Y*(u)  £  CmXm,  suc}i 


that 


M(ejujy 


X(u)  0 

0  I 


X{u)  0 
0  Y(  w) 


<  0 


holds  for  all  ui  £  [0,  27t],  and 

/  trace(y(m)) —  <  1 

Jo  2x 


(25) 

(26) 


This  condition  is  in  fact  very  similar  to  Condition  1  for  Robust  Hoo  performance.  The  only 
addition  is  the  incorporation  of  the  function  F (u),  which  can  be  assumed  continuous.  Heuristically, 
for  m  =  1,  F(w)  allows  for  the  gain  to  be  larger  than  1  at  some  frequencies,  provided  that  it  is 
compensated  at  other  frequencies  by  keeping  the  total  effect  fY(w)du>  less  than  1;  this  imposes  an 
“average  over  frequency”  performance  which  corresponds  to  the  H2  norm. 
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The  main  result  of  this  paper  is  that  this  test  answers  the  robust  71 2  performance  problem. 
As  in  the  case  of  Condition  1,  different  results  can  be  stated  in  accordance  with  the  nature  of  the 
perturbations  in  A,  which  exactly  parallel  Propositions  1,  2  and  3. 

For  the  first  one  involving  LTI  perturbations,  the  LFT  A*M  is  an  LTI  system  so  the  result  can 
be  stated  using  the  standard  definition  (13)  of  the  7L2  norm,  and  proved  with  elementary  frequency 
domain  tools. 

Theorem  6  Suppose  Condition  2  holds  for  matrix  functions  X(oj),Y(uj).  If  A  £  Balti,  then  the 
system  is  robustly  stable  and 

sup  ||A*M||2<1  (27) 

AeBAi/n 


Proof:  The  first  block  of  the  inequality  (25)  gives  ||X(o;)2  Mn(e3W)I(u)  2||oo  <  1,  which  implies 
(see  [26])  robust  stability  of  the  system  under  LTI  perturbations.  Furthermore,  defining 

X~i  0 

0  I 

we  conclude  that  for  some  e  >  0,  and  all  u>, 


M 


X  2  0 
0  I 


M 


I  0 


(29) 


Fix  A  G  Ba,  LTI.  For  any  fixed  frequency,  since  A(eiu),  X^(u)  commute,  we  can  replace  M  by 
M  in  Figure  1,  giving  A (e^u)-k  M(e^w)  =  A(eJW)  *  M(e^).  Using  (29),  we  have 


Is/HI2  +  k(^)|2  <  (1  -  c)|p(w) I2  +  u(u)*Y (oj)u(lu)  -  e|,u(a;)|2  (30) 


where  we  use  the  signal  denominations  of  Figure  1. 

Since  A  is  LTI,  contractive  we  have  \p(u)\2  <  \q(u>)\2:  which  leads  to 

u(lj)*(A  -k  M )(cu)*(A  *  M)(u>)u(u)  =  \y(u)\2  <  u(u))*Y (lo)u(lo)  -  e\u(u)\2  (31) 
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Since  this  holds  for  any  u(oj),  we  have 


(A  *  M)(u)*( A  *  M)(co)  <  Y(u)  -  el  (32) 

across  frequency.  Computing  the  trace  and  integrating  gives,  using  (26), 

f  2n  j 

||A*M||2<  /  trace(Y(u)) - me  <  1  (33) 

Jo  2?r 

□ 

We  have  obtained  a  convex  sufficient  condition  for  Robust  H2  performance  under  LTI  uncer¬ 
tainty;  we  now  show  that  it  has  the  same  necessity  properties  as  Condition  1  for  the  Hoo  case.  To 
state  the  following  results  for  which  include  non-LTI  systems,  we  adopt  the  approach  of  (24)  and 
give  the  following  definition: 

Definition  1  The  uncertain  system  ( M ,  A)  with  input  u  €  l™  has  robust  H2  performance  if  it  is 
robustly  stable,  and  there  exists  77  >  0  such  that 

sup  ||A*M||pym  <  1  (34) 

agba  " 


The  first  result  concerns  the  case  of  frequency- dependent  X-scales.  For  brevity  the  proof  is  not 
included  here,  although  we  will  later  remark  on  the  modifications  to  the  proof  of  Theorem  8  needed 
to  obtain  it.  A  full  account  is  given  in  [32]. 

Theorem  7  There  exists  v  >  0  such  that  the  system  (M,  A)  has  robust  H2  performance  for  A  e 
Ba"  if  and  only  if  there  exist  bounded  variation  functions  X (u)  G  X,  y(u>)  satisfying  Condition 

2. 

Theorem  7  gives  indication  that  there  is  mild  conservatism  involved  in  using  Condition  2  for 
LTI  uncertainty,  in  a  totally  analogous  way  to  Proposition  2  for  Ttoo  performance. 


14 


In  the  Hoo  case  there  was  also  supporting  empirical  evidence  with  computation  of  lower  bounds 
for  the  LTI  case  based  on  /i  [26],  which  is  not  available  for  'H-i-  In  fact,  the  restriction  on  causality  of 
the  LTI  perturbations  will  provide  an  additional  gap  for  H2  performance.  To  see  this,  consider  the 
case  of  unstructured  uncertainty,  where  the  TLoo  conditions  are  known  to  be  exact  [26].  In  the  7f2 
case,  with  scalar  inputs,  it  is  easy  to  show  that  Condition  2  is  exact  for  non-causal  LTI  perturbations, 
by  simply  choosing  A(eja')  to  produce  the  worst  gain  at  every  frequency;  this  interpolation  is  in 
general  only  possible  with  non-causal  A.  The  gap  due  to  causality  has  not  been  quantified  in 
general,  but  the  results  of  [45]  (see  Section  6)  suggest  that  it  is  not  significant. 

In  any  case,  the  only  necessary  conditions  available  for  the  Hoo  and  the  H2  frequency  dependent 
scales  tests  are  Proposition  2  and  Theorem  7,  both  indicating  that  these  gaps  are  a  modest  price 
to  pay  for  a  convex  characterization. 

We  now  turn  to  the  constant  X-scales  condition: 

Theorem  8  The  system,  (M,  A)  has  robust  7f2  performance  for  A  €  Baltv,  if  and  only  if  there 
exists  a  constant  matrix  X  E  X ,  and  a  bounded  variation  function  Y(u>),  satisfying  Condition  2. 

Proof:  For  simplicity,  the  proof  will  be  described  in  detail  for  the  case  of  scalar  inputs  u  £  l\, 
and  for  uncertainty  A  =  diag[Ai, . . . ,  Ara]  consisting  only  of  full  blocks.  For  the  general  case  see 
the  remarks  at  the  end  of  the  section. 

[Sufficiency]:  The  first  block  of  (25)  gives  \\X%  Mii(^“)X"2  ||oo  <  1,  which  implies  [26,  38]  robust 
stability  of  the  system  under  LTV  perturbations.  Also,  J2  and  A  commute,  so  define  M  as  in 
(28),  which  verifies  (29),  and  leads  to  (30),  which  can  be  integrated  across  frequency  to  give 

/27r  rh  1 

\\y\\2 +  \\q\\2  <(1- €)\\p\\2  +  u(u)*Y(u)u(u)  —  -  e||u||2  (35) 

Since  ||A||  <  1,  then  |[p||  <  ||g||,  leading  (for  scalar  u)  to 

J[2ir  1 

'  Y(u)\u(u)\2- - e|M|2  (36) 

0  2tt  '  ’ 
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(37) 


Fix  r/  >  0;  for  u  £  Wv,  we  invoke  Lemma  5  to  bound 

r 2tt  j  ,  r In  J.  . 

Jq  Y{u)\u{u)\^<VTV{Y)  +  \\utf  Y(u)  — 

Substituting  (37)  into  (36)  and  using  (26)  leads  for  ||«||  <  1,  to 

||(A  *  M)u\\2  <  tjTV(Y)  +  1  —  e  (38) 

By  choosing  small  enough  t ],  this  implies 

SUp  [|(A-*rlW)||M/I)  <  1 

A£Ba 

[Necessity]:  The  converse  implication  is  based  on  an  extended  “S-procedure  losslessness”  theorem 
on  quadratic  functions  on  1%.  This  type  of  result  was  first  obtained  by  Megretski  and  Treil  in  [25] 
for  a  finite  number  of  scalar  quadratic  forms.  The  nature  of  the  constraints  defining  the  white  noise 
sets  Wv  requires  the  extension  of  this  procedure  to  quadratic  functions  on  which  take  values  on 
the  function  space  Cr[ 0,27t]. 

Let  2  =  col(z!, . . . ,  Zn+i)  be  the  vector  of  all  inputs  to  the  M  system,  where  z\  . . .  zn  partition 
p  in  correspondence  with  the  blocks  Ai , . . . ,  An,  and  zn+1  =  u.  Analogously  ( Mz)i ,  i  =  1 . . .  n  +  1 
denotes  the  partition  of  the  output  of  M. 

Now  define  the  following  scalar  valued  quadratic  functions  of  z  £  I2, 

o-i(z)  =  \\(Mz)i\\2  -  ll^ll2,  i  =  l...n  +  l  (39) 

Roughly,  the  motivation  behind  these  functions  is  the  following:  if  the  crt  are  all  non-negative 
at  a  certain  2  /  0,  then  M  expands  this  signal  in  all  the  channels,  and  therefore  a  contractive, 
structured  LTV  operator  A  can  be  constructed  where  the  Hoo  performance  is  violated.  Now  for 
this  signal  z  to  violate  77 2  performance,  the  u  portion  must  be  white  (belong  to  Wn  as  defined 
in  Section  3);  this  requirement  is  imposed  by  an  additional  quadratic  constraint.  Consider  the 
function  p  :  I2  1-+  Cr{ 0,  2i r], 

p(z)  =  Fzn+ 1  ( 5 )  (40) 
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where  FZn+1  is  defined  as  in  (15).  With  this  definition,  zn+i  £  Wv  if  and  only  if  p(z)  £  f?(0, 77),  the 
ball  of  radius  77  in  £#[0,  2ir], 

By  considering  the  real  Banach  space  V  =  Rn+1  ©  Cj?[0,  27t],  we  can  collect  all  these  functions 
together  in  a  quadratic  map  A  :  1 2  V,  given  by 

A(*)  =*  (<ri(z), an+1(z),  p(z))  (41) 

Define  a  set  in  V, 

V  =  {A(*):*£/2,||*||  =  I}  (42) 

The  previous  discussion  suggests  that  robust  performance  will  be  violated  if  for  all  77  >  0, 
V  intersects  the  set  K,  :=  {(ri, . . . ,  7-n+1,  g)  :  r,  >  0,  ^(s)!^  <  77}.  For  technical  reasons,  we 
introduce  instead  the  set 

fCt:={(ru...,rn+1,g):  r;  >  -e2,  ||^(s)||^  <  77}  (43) 

We  formalize  these  ideas  in  the  following  Proposition,  which  reduces  robust  performance  to 
geometric  separation  condition  in  the  space  V. 

Proposition  9  Suppose  ( M ,  A)  has  robust  H2  performance  for  A  £  Baltv.  Define  V  as  in  (42), 
K.e  as  in  (43)  (for  77  given  in  Definition  1).  Then  there  exists  e  >  0  such  that  V  fl  /C£  =  0. 

To  bring  in  the  Hahn-Banach  Theorem,  we  note  that  K.t  is  open  and  convex  in  V,  and  that 

Proposition  10  The  closure  V  of  V  is  convex  in  V. 

Propositions  9  and  10  are  proved  in  the  Appendix.  By  choosing  K,\  =  V,  IC2  =  /C£,  we  are  in  a 
position  to  apply  Theorem  4,  and  obtain  the  corresponding  P  £  V*,  T  ^  0,  a  £  M.  The  structure 
of  V  and  the  Riesz  representation  theorem  imply  that  P  can  be  represented  by  (aq, . .  .,xn+i,  W), 
where  £  R,  $  £  BV[0,2tt}.  Then  (12)  yields 

/*27T 

y2xiai(z)+  /  p(z)dm  <  a ,  \/zel2,  ||*||  =  1  (44) 

i=i  Jo 
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n + 1  p  2w  (  2 

g>i+yo  9»>a,  (45) 

Concentrating  on  (45),  we  conclude  that  Xi  >  0,  i  =  1 . .  .n  +  1;  also,  since  /Ce  contains  a  ball  of 
0,  and  f  ^  0,  then  a  <  0.  Now  turning  to  (44),  it  is  possible  to  perturb  Xi  to  make  them  strictly 
positive,  and  since  <Ji(z )  are  bounded  functions,  (44)  will  still  hold  for  a  new  value  a  <  0.  Similarly, 
xn+\  can  be  normalized  to  1. 

It  only  remains  to  rewrite  (44)  using  the  definitions  of  ct;,  p.  In  the  first  place,  simple  manipu¬ 
lations  show  that 


where  X  =  diagfzj/, . . .  ,xnI ]  >  0  is  a  scaling  which  commutes  with  A.  Secondly,  an  integration 


by  parts  (note  FZn+1  (0)  =  FZn+1  (2? r)  =  0)  gives 

f2n  r2v  t-2  ir  j 

l  FZn+i{s)<M{s)  =  -  jf  *(s)F^+1(s)ds=-  jo  ^)(|zn+1(S)|2-||^+1||2)-  (47) 

Defining  Y (s)  =  1  +  \P(s)  -  J02?r  <H(,s)|^,  the  right  hand  side  of  (47)  is  equal  to 

(yW-1)K«WI2^  =  ([o  i-y\z’z)  <48) 

Combining  everything  into  (44),  we  obtain 

((V  *  °r  M-  *  °  Vzel2,\\z\\  =  l  (49) 

This  implies  (25)  holds.  Finally,  from  the  definition  of  Y(s),  we  know  that  JQ27r  Y(s)^  =  1,  and  a 
small  perturbation  in  Y  will  preserve  (25)  and  yield 

f2vr  ds 

L  Y(e)2i  <  1  (50> 

which  is  (26)  for  this  scalar  case. 

□ 
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We  now  comment  briefly  on  the  various  extensions  to  the  above  proof;  these  are  developed  in 
detail  in  [32]. 

Multivariable  noise. 

If  u  £  Z™,  then  from  (23)  the  performance  quadratic  constraint  is  changed  to 

an+1(z)=\\(Mz)n+1\\2-^\\zn+1\\2  (51) 


and  (22)  indicates  the  natural  definition  for  p(z)  =  FZn+1(s),  given  by 


fS  n 

F*n+1(S)  =  /  Zn+1(u)zn+1(u)*—  -  —  l^n+lll^m 

J  o 

FZn+1(s)  now  takes  values  in  the  space  of  continuous,  hermitian  matrix- valued  functions;  the  dual 
of  this  space  can  be  identified  with  the  space  of  hermitian,  bounded  variation  matrix  functions 
on  [0,  27t]  (up  to  a  constant  matrix),  with  the  convention 


U(g)  = 


trac  e(g(s)d'9(s)) 


The  proof  then  follows  in  a  similar  way,  giving  F(s)  =  \P(.s)  +  ^1  -  /Q27r  traceU/(s)|^  F  which 
satisfies  (25)  and  is  perturbed  to  satisfy  (26). 


61  perturbations  in  A. 

If  the  i-th  block  of  A  is  6Iri,  then  the  scalar  quadratic  function  cr,-  must  be  replaced  (see  [32]) 
by  a  matrix- valued  function 

Si(*)  =  JQ  l{Mz)i{u)(Mz)*i{u)  -  *(«>*(«)]  ^ 

which  takes  values  in  the  space  of  hermitian  T{  x  r,  matrices.  The  functionals  in  this  space  are  of 
the  form  Tj^ (A)  =  trace(XjA),  where  Xt  is  a  full,  hermitian  matrix.  The  argument  then  proceeds 
in  a  similar  fashion,  X*  becoming  a  sub-block  of  the  scaling  matrix  X. 
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Slowly  varying  perturbations 

To  modify  this  proof  to  obtain  Theorem  7,  we  need  a  quadratic  constraint  for  slowly  varying 
perturbations.  It  can  be  shown  (see  [35,  30,  32])  that  given  p,  q  G  I2,  if 

fs+h  rlhJ 

/  (k!2-M2)^>o 

for  every  interval  of  length  h,  then  there  exists  a  contractive  operator  A,  with  ||AA  —  AA||  <  v  = 
2sin(|),  such  that  A q  =  p.  Hence  it  is  natural  to  replace  the  scalar  (7{  with  a  quadratic  function 
Vi  ■  h  Cfi[0,27r], 

ps+h  1 

[<w(*)](*)  =  js  (\(Mz)i\2  —  \zi\2)^ 

The  corresponding  functional  will  be  a  bounded  variation  function  Xi(s),  which  will  give  the  fre¬ 
quency  varying  portion  of  the  X  scale.  In  this  way  one  can  also  handle  a  combination  of  slowly 
varying  and  LTV  perturbations  for  robust  H2  performance,  as  was  done  in  [30]  for  'Hoc  performance. 

5  Computational  Issues 

A  test  has  been  developed  in  Section  4  which  characterizes  robust  7^2  performance  analysis  of  an 
uncertain  system.  This  test  is  an  infinite  dimensional  convex  feasibility  condition  on  the  unknowns 
X  and  Y,  specified  as  a  Linear  Matrix  Inequality  (LMI)  across  the  frequency  axis,  of  a  similar 
complexity  as  Condition  1  for  robust  TL^  performance. 

There  are  two  standard  approaches  for  handling  the  infinite  dimensionality  of  these  conditions, 
and  turn  them  into  finite  dimensional  LMIs,  for  which  efficient  algorithms  are  available  [6]:  one 
used  in  [4]  is  to  grid  the  frequency  axis,  the  other  is  to  select  a  finite  set  of  basis  functions  and 
search  for  a  scaling  in  the  span  of  these  functions,  which  reduces  to  a  single  LMI  via  a  state-space 
approach  (see,  e.g.,  [3]). 

Both  approaches  can  indeed  be  applied  to  Condition  2,  and  involve  minor  modifications  to  their 
counterparts  for  Condition  1.  We  demonstrate  this  by  commenting  on  the  gridding  approach  for 
this  problem:  Condition  2  is  approximated  by  considering  frequency  points  0  =  Wq  . .  .wn  =  27r. 
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Although  this  approximation  offers  no  hard  guarantees,  since  it  is  based  on  the  frequency  domain  it 
allows  for  engineering  judgement  to  be  used  in  choosing  the  number  and  location  of  the  grid-points. 
The  finite  dimensional  approximation  to  Condition  2  is  the  LMI  problem 

M(^r  *  J  M(ui)  -  X0l  °  <0  i  =  l...N  (52) 

Xi  >  0  i  =  l...N  (53) 

N 

£  trace(K')(^i  ~  1)  <  1  (54) 

2  =  1 

where  the  unknowns  Y{  are  hermitian  matrices  and  the  Xt  structured  matrices.  For  the  LTV  test, 
=  X  is  constant  across  the  ojt ,  which  makes  conditions  (52-54)  intrinsically  coupled  across 
frequency.  For  the  LTI/slowly  varying  test,  we  use  different  variables  X, ,  i  —  1 . .  .N.  Although 
(54)  still  involves  all  frequency  points,  the  following  strategy  can  be  used  to  decouple  the  problem 
across  frequency: 

•  For  each  fixed  frequency  point  u>i,  pose  the  problem: 

Minimize  trace(lj)  ,  subject  to 

[*‘  "]  <0  (55) 

Xi  >  0  (56) 

The  problem  of  minimizing  a  linear  function  of  the  unknowns,  subject  to  an  LMI  constraint, 

falls  in  the  class  of  eigenvalue  problems  (EVPs)  considered  in  [6],  and  can  be  computed 
efficiently. 

•  Given  all  the  solutions  Fj  . .  .Yn,  compute  ]TL_i  trace(lj)(o;j- —  u>i-i),  and  compare  the  answer 
to  1.  More  directly,  this  sum  will  provide  an  approximation  to  the  square  of  the  worst-case 
Ti2  norm  of  the  system;  this  follows  from  the  fact  that  to  test  if  the  worst-case  'H2  norm  is 
less  than  7,  it  suffices  to  change  1  for  72  in  (26)  or  (54). 
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6  Connections  to  Mixed  ^2/^00  Performance 


In  this  section  we  relate  these  results  to  earlier  work  in  the  so-called  mixed  ^2/^00  problem.  There 
are  many  versions  of  this  problem  in  the  literature  (a  few  are  [5,  22,  45,  34]),  all  of  which  attempt  to 
get  a  handle  on  robust  77  2  performance  by  studying  first  the  situation  where  there  is  no  uncertainty, 
but  the  performance  specification  is  a  combination  of  the  772  and  Hoo  norms. 

A  mixed  'H2/W00  performance  problem  can  in  fact  be  cast  naturally  in  our  setting,  and  leads 
to  an  analysis  test  which  looks  exactly  like  Condition  2,  except  that  the  scaling  matrix  X(ov)  does 
not  appear  and  is  fixed  to  be  the  identity. 


Proposition  11  Consider  a  system  M  —  [M\  M2]  where  the  input  z  is  partitioned  in  the  vectors 
z\,  Z2  E  Cm.  The  following  are  equivalent: 


0): 

(ii)  : 


3p  >  0  :  sup  ^  \\Mz\Y  :  \\zt\\z  +  —  \\z2\\z  <  1,  e  Wn  f  <  1 

TYl 


There  exists  Y(u)  =  Y(co)*  : 
M(ejuj)*M(ej“)  - 


I  0 

0  Y (u) 


<0  V  u  £  [0, 27r] 


f2n  /T  rt  \\du) 

J ^  trace(F(o;))—  <  1 


(57) 


(58) 

(59) 


This  result  is  proved  along  the  same  lines  as  Theorem  8,  only  that  since  there  is  no  uncertainty, 
one  only  considers  the  whiteness  constraint  p(z)  =  FZ2 ,  and  the  performance  quadratic  constraint 
a(z)  =  ||A/A||2  —  ||^i||2  —  ^ 1 1 1 1 2 >  with  ||^2||2  weighted  by  ~  for  the  reasons  explained  in  Section  3. 

Conditions  (58)  and  (59)  are  therefore  interpreted  as  a  mixed  performance  problem  where  a 
portion  of  the  input  signal  is  constrained  to  be  white.  Various  problems  like  this  were  considered 
in  [45]  for  continuous  time,  with  different  assumptions  on  the  relationship  between  z\  and  z2.  It 
is  shown  in  [45]  that  the  performance  costs  are  not  substantially  different  for  these  alternatives, 
and  subsequently  the  attention  is  concentrated  on  the  case  where  z\  is  restricted  to  be  causally 
dependent  on  z^-  State-space  methods  for  both  analysis  and  synthesis  for  this  alternative  are  given 
in  [45]  and  the  sequel  paper  [16]  (a  stochastic  version  appears  in  [34]). 
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Our  condition,  in  contrast,  corresponds  to  the  case  where  there  is  no  such  causality  restriction, 
which  is  only  treated  summarily  in  [45].  While  at  the  level  of  the  mixed  performance  problem  it  is 
not  obvious  which  alternative  to  prefer,  the  version  considered  here  has  advantages  from  the  point 
of  the  robust  'Hi  problem,  since  it  is  allows  for  the  inclusion  of  X  scales  in  a  convex  condition  with 
the  strong  interpretation  given  in  Theorems  7  and  8. 

Remarks: 

•  These  causality  issues  are  strongly  connected  to  the  remarks  of  Section  4  regarding  the  causal¬ 
ity  of  LTI  perturbations.  A  possible  conclusion  from  the  results  in  [45]  is  that  gaps  due  to 
this  are  not  very  significant,  although  this  issue  warrants  further  investigation. 

•  Analogously  to  the  case  treated  in  detail  in  [45],  the  mixed  conditions  (58-59)  reduce  to  a 
finite  dimensional  test  if  M(e^)  is  rational.  In  fact,  a  Schur  complement  operation  and  some 
algebra  shows  that  (58)  is  equivalent  to 

HAfiHoo  <  .1  (60) 

M%{I  ~  M1M$)~1M2  <  Y  Vome[0,2tt]  (61) 

This  implies  that  (58-59)  can  be  tested  by  first  checking  (60),  and  then  imposing  that 
||A-1M2||2  <  1,  where  N  is  the  inversely  stable  spectral  factor  (see,  e.g.,  [46])  satisfying 
I  —  M\M j  =  NN*.  Both  these  operations  can  be  computed  efficiently  by  the  same  state- 
space  techniques  used  in  the  mixed  Hi/Hoo  literature. 

7  Robust  7Y2  Synthesis 

Having  obtained  conditions  for  robust  Hi  performance  analysis  under  structured  uncertainty,  it  is 
natural  to  consider  the  problem  of  controller  synthesis.  If  the  nominal  system  M  is  obtained  as 
the  closed  loop  in  a  feedback  configuration,  the  problem  is  to  design  the  controller  such  that  M 
satisfies  the  robust  Hi  performance  conditions. 
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It  is  unlikely  that  a  tractable  global  solution  will  appear  for  this  problem,  since  none  is  known 
for  the  case  of  performance.  Except  for  very  special  configurations,  the  only  general  method 
for  robust  7ioo  synthesis  is  the  so-called  “D-K”  iteration,  where  an  an  analysis  step  (Condition  1) 
is  alternated  with  Hoo  synthesis. 

Such  iteration  schemes  can  easily  be  extended  to  robust  H2  synthesis,  as  is  now  described. 
Assume  that  the  functions  X(u;)  and  Y{ui )  are  rational  and  satisfy  Condition  2  (they  could  be 
obtained  by  fitting  frequency  points,  or  with  basis  functions).  They  are  both  positive  definite 
across  frequency,  so  by  a  spectral  factorization  (see  [46])  they  can  be  expressed  as 

X(u)  =  D(ejw)*D(ejw),  Y(u)  =  E(e^)*E{e^),  (62) 


where  D,  E  are  rational,  stable,  stably  invertible  transfer  functions,  and  D  has  the  structure  of  X. 
Using  these  factorizations,  it  follows  that  (25)  can  be  rewritten  as 


•  For  fixed  D,  E,  reduce  the  norm  in  (63)  by  Ttoo  synthesis. 


•  For  a  fixed  controller,  solve  the  analysis  problem  for  D,  E. 


As  for  the  standard  D-K  iteration,  each  step  in  this  “(D,E)-K”  iteration  can  be  shown  to  improve 
the  robust  performance  cost,  but  there  is  no  reason  to  expect  convergence  to  the  global  optimum. 

The  previous  iteration  was  based  in  H0 o  synthesis.  An  alternative  is  suggested  by  the  discussion 
in  Section  6,  where  at  the  synthesis  step  one  only  includes  the  D  scales  with  the  plant,  and  employs 
the  design  schemes  for  the  mixed  H2IH00  problem.  As  remarked  before,  the  techniques  in  [16] 
correspond  to  a  slightly  different  mixed  problem,  and  it  is  not  clear  whether  they  extend  to  the 
situation  of  Proposition  11.  The  approximate  method  presented  in  [8]  could  also  be  used. 

Therefore  a  number  of  issues  remain  open  for  future  research,  regarding  this  second  iteration 
and  the  comparison  between  the  two  alternatives  for  practical  problems. 
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8  Conclusion 


The  results  in  this  paper  restore  the  7f2  performance  paradigm  to  the  level  of  mainstream  robust 
control:  there  is  no  longer  a  significant  advantage  from  the  point  of  view  of  robustness  analysis 
in  the  consideration  of  Tioo  norms,  since  it  is  essentially  no  harder  to  test  for  the  corresponding 
^-performance  problem. 

Condition  2  appears  in  fact  as  a  summary  of  tractable  exact  conditions  for  robustness  anal¬ 
ysis.  Setting  the  blocks  in  X  to  be  either  constant  or  frequency-varying  selects  between  LTV 
or  LTI  (slowly-varying)  uncertainty,  or  a  combination  thereof.  Selecting  either  a  constant  or  fre¬ 
quency  varying  Y  chooses  between  maximum  over  frequency  or  average  over  frequency  performance 
(roughly,  'Hoc  or  H2  performance;  in  rigor,  for  Tioo  the  trace  must  be  removed  from  (26)).  Also, 
combinations  of  H2  and  Hoc  performance  can  be  studied  in  the  style  of  Section  6,  by  including  Y 
terms  only  for  the  signals  which  are  assumed  white.  For  all  these  choices,  we  have  an  exact  char¬ 
acterization  of  the  robust  performance  problem  for  which  Condition  2  is  necessary  and  sufficient. 

A  number  of  research  questions  are  raised  by  these  results.  One  of  these  is  the  issue  of  com¬ 
putation  of  Condition  2;  we  demonstrated  its  tractability  in  Section  5,  but  further  research  and 
practical  experience  with  algorithms  are  in  order.  Another  important  question  which  we  will  in¬ 
vestigate  in  future  papers  is  the  applicability  of  state-space  methods  for  this  problem,  which  would 
reduce  Condition  2  to  a  finite  dimensional  LMI  in  the  case  of  constant  scales.  Finally,  the  problem 
of  controller  synthesis  for  robust  7f2  performance  is  now  reduced  to  the  complexity  level  of  robust 
Hoc  synthesis;  further  investigation  is  required,  in  particular  on  the  different  iteration  schemes 
which  where  proposed  in  Section  7. 

To  some  extent,  this  paper  closes  a  cycle  of  research  which  originated  in  the  70s.  The  methods 
developed  during  this  period  have  reached  the  maturity  to  address  one  the  main  problems  which 
motivated  the  appearance  of  robust  control  theory.  The  impact  of  this  theory  goes,  however,  far 
beyond  the  robustness  of  LQG  regulators,  constituting  a  fundamental  addition  to  the  understanding 
of  models,  uncertainty,  and  feedback  control. 
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Appendix 

We  will  give  proofs  of  Propositions  9  and  10. 

Lemma  12  Let  A  :  l2~*V  be  defined  as  in  (41).  If  z,  f  £  I2,  and  X  is  the  delay  operator,  then 

\\z  +  ^kf\\2  k—->  IMP  +  ||/||2  (64) 

A(z  +  Xkf)k^?A(z)  +  A(f)  (65) 

where  (65)  means  convergence  in  the  topology  ofV  —  Eri+1  0  Cr[0,2tt], 


Proof:  The  main  observation  (which  is  easily  seen  by  looking  at  the  shift  operator  in  the  time 
domain)  is  that  ( z ,  A k f)  0  for  any  functions  z,f£  /2.  From  this  it  follows  that 

\\zi  +  XkM2  k H^ll2 +  ||/4'||2  (66) 

I HMzfi  +  (MXkf)i\\2  ^  \\(M z)i\\2  +  UMffiW2  (67) 


where  in  (67)  we  use  MXk  =  A kM  from  the  time  invariance  of  M.  This  implies  (64),  and  also 


ofiz  +  A  kf)  k=^'  Oi(z)  +  off)  i  -  1 . .  .n  +  1 


(68) 


We  now  show  that 


p(z  +  A kf)  p(z)  +  p(f) 


(69) 


with  convergence  in  the  sense  of  Cr[ 0,2tt].  Starting  from  (40)  and  (15),  some  algebra  gives 

.  du  s 


(70) 


p(z+Xkf)  -  p(z)  -  p(f)]  (s)  =  (jf  2  Re[z*n+1(u)  ej“k  fn+1(u)]~  -  S-Re  (zn+1 ,  Xk  fn+1) 

It  therefore  suffices  to  show  that  the  sequence  of  functions 

Ms)  :=  jf  zn+ i(w)  eJU,fc/„+ i(w)^  (71) 
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converges  uniformly  to  0  as  k^oo.  Pointwise  convergence  follows  from  tpk(s)  =  {l[o,s]^n+i>  Afc/n+i)- 
If  convergence  were  not  uniform,  we  could  find  e  >  0,  a  subsequence  kj  and  points  skj  with 
\<Pkj(skj)\  >  £■  By  compactness,  taking  a  partial  subsequence  we  can  assume  skj  'J— ?  -So-  Now 

<Pkj(sk3)  =  <Pkj(s0)  +  J'  3  «*+i(w)  e>wk*  fn+1(v)^  (72) 

so 

0  <  e  <  \<Pk3(skj)\  <  |^(s0)|  + 

The  right  hand  side  of  (73)  converges  to  0  from  the  pointwise  convergence  of  yu. ,  and  the  fact  that 
|2n+i||/n+i|  £  £i[0,27t].  This  is  a  contradiction,  so  we  have  shown  (69),  which  together  with  (68) 
implies  that  A(z  +  A kf)  *-72?  A(z)  +  A (/)  in  the  topology  of  V. 


! 
J  Si 


Sk,  dto 

\zn+i\\fn+i\  — 


(73) 


□ 


Remark:  The  previous  Lemma  can  easily  be  extended  to  the  sum  of  shifted  versions  of  N  signals 
z(°\  . .  G  h,  giving 


JV-l  N- 1 

ii  \krz^f y  ik(r)n2 

r= 0  r= 0 

/N- 1  \  N-l 

a  (  y  xkTz(r)  A^(r)) 


v  r=0 


r= 0 


(74) 

(75) 


Proof  of  Proposition  10 


Let  co(V)  denote  the  convex  hull  of  V;  an  element  Ao  E  co(V)  is  a  convex  combination  of  the  form 

N-l  N-l 

Ao=Y  arHz(r) ),  «r  >  0,  J]  =  1,  ||^(r)||  =  1  (76) 

r=0  r=0 

Define  fk  =  ^/a^Xkr  z^rh  From  (75)  and  the  quadratic  nature  of  A,  it  follows  that 

N-l 

A(fk)  ^  2  A(V^W)  =  Ao  (77) 

r=0 
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From  (74),  we  find  ||/fc||2  * ar  =  1.  Therefore  A( jjj^yy)  A0,  so  A0  E  V. 
We  have  shown  co(V)  C  V.  This  implies  that  co(V)  C  co(V)  C  V,  so  V  is  convex. 


□ 


For  the  next  Lemma  we  consider  the  configuration  of  Figure  2,  where  a  disturbance  signal  d  is 
injected  at  the  output  of  A. 


Figure  2:  Uncertain  system  with  injected  d 


Lemma  13  Given  e  >  0,  suppose  V  fl  Kt  ^  0.  Then  there  exists  an  operator  A  G  Baltv,  and  I2 
signals  z  =  col(p,u),  Mz  =  col(q,y),  d,  satisfying  the  equations  of  Figure  2,  with  u  G  Wv,  ||z||  =  1, 
||d||  =  0(e),  and  ||y||  >  ||«||  -  0(e). 

Proof:  By  hypothesis,  we  can  find  £  =  col(v,u),  ||^||  =  1,  such  that  A(^)  E  ICe.  This  means 
that  p(z)  G  B( O,??),  so  u  E  Wv,  and 

(Ti(z)  =  \\(Mz)i\\2  -  ||^j|2  >  -e2,  i  =  1 . . . n  +  1  (78) 

From  (78)  we  can  find  a  contractive  operator  A*  :  and  a  disturbance  di,  ||di||  =  0(e),  such 

that 

Z{  =  A i(Mz)i  +  di  i  =  1 . . .  n  (79) 

Also,  ||y||  >  ||n||  —  e  follows  from  the  constraint  on  an+i-  This  would  complete  the  proof,  except 
that  the  given  operator  A  =  diag[Ai, . . . ,  A„]  need  not  be  causal.  For  (79)  to  hold  with  causal  A 4- 
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would  require  (see  [36])  the  stronger  condition 


\\PT (M z)i\[2  -  ||PT^||2  >  -0(e)  VT  G  Z+  (80) 

(without  loss  of  generality,  assume  2  supported  in  Z+)  where  PT  denotes  the  truncation  operator: 

for  i  >  T  <81) 

In  fact,  (80)  can  indeed  be  satisfied  by  means  of  a  construction  due  to  Shamma  [38],  where  z  is 
replaced  by  a  signal  z  obtained  from  repetition  of  shifted  versions  of  z.  Consider  the  signal 


JV-l 


(82) 


T— 0 


Choosing  =  0(e),  and  sufficiently  large  k,  it  can  be  shown  that  this  signal  satisfies  (80)  as 
required.  The  details  of  this  argument  are  quite  involved  and  are  omitted  (refer  to  [32]  for  a  complete 
treatment  of  causal  perturbations).  Note  also  that  p{zk)  k^-$'  p(z )  from  (75),  and  \\zk\\  k--2?  1  from 
(74),  so  zk  will  satisfy  all  the  required  conditions  for  large  k. 


□ 


Proof  of  Proposition  9 

From  robust  stability,  (/  —  AMn)'1  exists  for  each  A  G  Baltv.  We  will  use  the  uniform  bound 

sup  ||(J-AM11)-1||  =  /3<oo  (83) 

AGBaltv 

which  appears  to  be  slightly  stronger  but  can  be  shown  (see  [32])  to  be  equivalent  to  robust  stability. 
We  also  know  by  hypothesis  that 

sup  ||A*M|U  =  7  <  1.  (84) 

AgBaltv 

Given  e  >  0,  suppose  K,e flV  ^  0.  We  apply  Lemma  13  and  construct  the  corresponding  A  G  Baltv, 
z  —  col(p.  u),  and  d.  We  now  state  the  following  bounds: 

-0{e)  +  ||u||  <  ||ff||  <  /3||M2i||||t/|j2  +  Tll^lla  (85) 
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The  lower  bound  is  a  direct  consequence  of  Lemma  13.  The  upper  bound  is  obtained  by  writing 
y  as  the  superposition  of  the  contributions  of  the  inputs  d  and  u  £  Wv  in  Figure  2,  and  using  the 
bounds  (83)  and  (84). 

If  fCc  fl  V  ^  0  holds  for  every  e,  we  can  choose  sequences  Ak,  zk  =  col (pk,uk),  ||^fc||  =  1, 
dk  corresponding  to  0.  From  the  bounds  (85)  and  7  <  1  it  follows  that  we  must  have 

||wfc||  kZl??  0;  also,  ||dfc||  =  O(ek)  0.  Now  from  (83)  the  gain  from  (d,u)  to  2  can  be  uniformly 
bounded  across  A,  which  implies  that  \\zk\\  k— ?  0,  which  is  a  contradiction.  Therefore  there  must 
exist  e  >  0  such  that  K,t  fl  V  =  0. 

□ 
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